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Abstract 



' (— I ' ' By using a formal analogy between statistical mechanics of mean field spin systems and 

Q^l analytical mechanics of viscous liquids -at first pointed out by Francesco Guerra, then re- 

cently developed by the authors- we give the thermodynamic limit of the free energy and 
the critical behavior of Curie Weiss models for a certain class of generaHzed spin variables. 
' Then, with the same techniques, we give a complete picture of the bipartite Curie- Weiss 

model, dealing with the same class of generaHzed spins. Ultimately we analyze further the 
existence of a minmax principle for the latter which mirrors the standard variational principle 
of canonical thermodynamics when generalized to multiple interacting parties. 

■ Keywords: MSC35, MSC76, MSC82 

'nJ- : 1 Introduction 

O \ The investigation of statistical mechanics of mean field spin systems is experiencing an increas- 

ing interest in the last decades. The motivations are two-fold: from one side, at the rigorous 
mathematical level, a clear picture is still to be achieved (it is enough to think at the whole 
community dealing with the case of random interactions as in glasses [E]), at the applied level, 
^ ' these toy models are starting to be used in several different context, ranging from quantitative 

■ sociology ^ to theoretical immunology [13j . 
It is then obvious the need for always stronger and simpler methods to analyze the enormous 
amount of "variations on theme", the theme being the standard and simplest dichotomic Curie- 
Weiss model (CW) [I]. 

Recently, inspired by the pioneering work of Francesco Guerra [TJ, we paved a clear way to man- 
age models with self- averaging order parameters by using the CW prototype as a guide [6]. Here, 
at first, we apply our scheme to work out the single-party CW with general spins (i.e. continuous 
spins with compact support and symmetric probability measure t.9|) and we solve in all details 
its thermodynamics. Then we switch to the case of bipartite systems with generalized spins |5j, 
both offering a clear picture of the thermodynamics as well as a digression on the connection 
of the coupled self-consistent equations for these models with the existence of an underlying 
minmax principle [8]. 

Overall we covered both the ways of investigation: from one side our analysis is mathematically 
clear (no powerful but not fully rigorous methods as replica trick or saddle points are used). 
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1 



from the other side it is automatically ready for being implemented into applied scenarios, i.e. 
the generalized bipartite system extends competitions in decision making [TT] among two com- 
munities by allowing their constituent to assume "softer" viewpoints among each other (with 
respect to the accept/reject perspective), being the spins ruling their will continuous instead of 
dichotomic. 



2 The ferromagnet with continuously distributed spin variables 

We deal with a system made by N i.i.d. spin r.v. Uj, i = 1...N, with probability measure /x(cr), 
having the following properties: 



i) fi{(7) = fi{-a), 

i.e. it is symmetric; 



ii) 3L : V6 > Oj;:l^^'df,{a)f{a) = jZ^JiJiia) f {a) = 0, 

i.e. it has compact support [— L/2,L/2]. 

In particular, denoting with E^^ = dp,{ai)...dji{ajq) the expectation values with respect to 
the spin variables, we notice that from i) it follows that Eo-[cr] = 0, and from ii) that, for a 
given bounded function of spin /(cr), it has to be Eo-[/(o")] < -^^[supo-g[_£,/2,L/2] /('^)]- 

The spins interact each other, in the way described by the Hamiltonian H^{a,h) 

^ N,N N 
HN{cr,h) = -—^GiGj -h^ai. (1) 

Partition function, pressure and free energy per site are defined as usual as 

AN{l3,h) = llogZjv(/3,/i), 

Of course we are interested in calculating the value of the free energy in thermodynamic limit, 
i.e. for A'^ — > oo, for describing the thermodynamics of the model. We can also define Boltzmann 
states of our system for a generic function of the spins ^at, as 

/ ( E<^j^C/Ar((T)exp(-/3i?Ar(c7,/l)) 

^^^^"^^ = ^Mh) ■ 

It is useful to define also the following quantities: 

1 ^ 

^N = ^^(Th (3) 



N 

N 

N 
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respectively the magnetization of the system, and the self overlap of spin variables (of course we 
have trivially oat = 1 VA^ in the case of dichotomic spin). We always have trivially {un) < L'^- 
We can even express the Hamiltonian ([1]) in terms of (3) and It is 

Hn{(t, h) = -N{^mlf + hruN) + ^oat. 
This is the starting point of the next section. 



2.1 The free energy in the thermodynamic limit 

We will follow the approach described in [7] [6]. With this purpose, let us introduce the function 
ipN{x, t) : M X M+ ^ M defined as 



1 

'n 



log E^^ exp 



tN 



xNtun I . 



(5) 



This function plays the role of the pressure (up to a sign) of a model defined with a Boltz- 
mannfaktor exp + xA^mTv) depending by the two other parameters (x,t) (it is just 

a matter of names, since both t and /? as well as x and h have the same range of definition). 
Anyhow it is strictly related to the pressure of our model, as states the following 



Lemma 1. It is, uniformly in N 

\ipN{x,t) + AN{x,t)\ < O 



(6) 



Furthermore we have that ipiy{0,x) = Aj^{0,x) is Lipschitz-continuous \/N. 
Proof. Due to convexity of logarithm we get 



\ipN{x,t) + AN{x,t)\ = 



log 



E^^ exp {^m% + xNtun) e 
E^jy exp (^m^ + xNtun) 



^ 1 1 ^ t_a^ ^ LH 

< — logEo-e 2 < 

and ^ is proven. Furthermore, trivially it is ipNiO,x) = —An{0,x) = — logEo-[e^'^]. Again by 
a convexity argument, joint with the compactness of the support of the cr-distribution, it is 



\ipNix,0) - (Pn{xo,^ 



log 
log 



E,e^ 



E^e^o'^ 

]g g(a;-xo)o-gXoc 



< logE^el'^-^ol'' < L\ 



X — Xq\ 



hence ipiy{x,0) is Lipschitz-continuous. □ 

Remark 1. We have that in the thermodynamic limit A{(3, h) 



(7) 
□ 



-^{x = h,t = 13). 



Remark 2. From (B^ it is easily seen that in thermodynamic limit the definitions of state 
and the one built with Boltzmannfaktor ex.p (tNm'j^/2 + xNtu^^ do coincide (of course replacing 
{(3, h) with {t, x) ), thus we actually won't distinguish them in the following. 
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The main idea, for solving the thermodynamics encoded into our Hamiltonian, is to relate 
the statistical mechanics system to an effective mechanical one, in which we naturally identify 
X with space coordinate, t with time, and the function (pN{x,t) with the mechanical action. In 
order to do this we notice that ipN{x,t) satisfies the differential problem ^ 

dt(pN{x,t) + ^{d:^ipN{x,t) f - jj^d^ipN{x,t) = in M X (0, +CX)) , 
ip N {x, 0) = -logE^e'"' onMx{t = 0}. ^' 

This is a Hamilton-Jacobi equation with a vanishing dissipative term in the thermodynamic 
limit. Defined uj\/{x,t) = dx^PN{x,t) the velocity field, we notice that it corresponds to magne- 
tization of the finite size system in our parallelism [6]. We have that UN{x,t) satisfies a Burger's 
equation again with a mollifier dissipative term: 



f dtUN{x,t) +UN{x,t)dxUN{x,t) - jj^dluN{x,t) =0 inMx(0, +00) . . 

\ Uiv(x,0) = -E^cie^'^/E.^e^'^ on M x {t = 0}. ^ ' 

Thus the problem of the existence and uniqueness of the thermodynamic limit is here translated 
into the convergence of the viscous mechanical problem to the free one. We can use a theorem, 
that resumes a number of results obtained by Peter Lax [10] and assures the existence of the 
solution for free problem: 

Theorem 1. For a general differential problem 

dMx,t) + ^{dMx,t)f = mMx(0,+oo) . . 

ip{x,0) = h{x) onRx{t = 0}, ^ ' 

and 

dtu{x,t) + u{x,t)dxu{x,t) = m M X (0, +00) , . 

u{x,0)=g{x) onRx{t = 0}, ^ ^ 

where h{x) is Lipschitz-continous, and g{x) = h'{x) G , it does exist and it is unique the 
function y{x, t) : M x M such that 

t f X — I t fx — y{x,t)^'^ 



(^(rE,t)=min|-(^^J + h{y) ^ = - ^ j +h{y{x,t)) (12) 

is the unique weak solution of Iil0\} . and 

x-y{x,t) 

u{x,t) = (13) 

is the unique weak solution of Furthermore, the function x y{x,t) is not- decreasing. 

It is easily seen that Lax's theorem gives us the solution for the free energy of the model. In 
fact, if we put u{x,t) = —M{x,t) in the solution of the free Burgers' equation, and use F as a 
short label standing for "free", we get that the minimizing function is yix,t) = x + tM{x,t), and 
the action of the mechanical model reads off as 

^2 

ipF(.x,t) = -A{x,t) = —M'^{x,t) -logE^[exp((T(x + tM(x,t)))]. 
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Therefore, we have the following expression for the free energy per site of our models, 

1 

= ^M\h,P)-^logE,[exp{aih + f]M{h,m]- (14) 

Finally we must prove convergence of the viscous problem to the free one. To this purpose, 
we can state the following 

Theorem 2. The function 

ipN{x,t) = -^log^ J -^exp[-iV((x-y)V2t-logE,exp(Ty)] (15) 
does solve equation ^ and it is 

\(pN{x,t) - fF{x,t)\ < O ^-^ 

Furthermore, the function 

/ exp [-N {{x - yf /2t - logE, expay)] 

UN(x,t) = J (16) 

I^e^P[-N{{x-y)y2t-logE^e^pay)] 

does solve equation ^ and it is 

1 



\uNix,t) + M{x,t)\ < O 

The proof is exactly analogue to the one given in |6j, and ultimately due to the uniform 
convexity of the exponent in (jlSp and (|16p . that we have here by construction, so we will not 
report it here. 

Therefore we have proven the existence of the thermodynamic limit for free energy and 
magnetization of our model. 

As symmetry breaking are fundamental even in statistical mechanics, we want to report hereafter 
some other considerations about the existence and the properties of a phase transition in our 
analogy. 

2.2 Phase transition and shock waves 

In this section we deeply study properties of the free Burgers' equation for the velocity field ([TT]) 
(that we remind is the analogue of the magnetization). We can write the straight line trajectories 
of the free system (i.e. the system in thermodynamic limit): 

^ " ' (17) 

As usual in these cases, we can find a solution for the magnetization along characteristics [4]. 
It is 

- nix, t) = Mix, t) = IEj^exp [a(x + tM(x,t))] _ 

^ ' ^ ^ ' ' E^exp[a{x + tM{x,t))] ^ ' 
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Remark 3. Putting {x = h,t = (3) in hi 9^) we recover the generalized self consistence equation for 
the magnetization. In particular, by choosing fj,{a) = {l/2)[6{a + 1) +S{a — 1)], we immediately 
recognize the well known hyperbolic tangent of the dichotomic CW model. 

An important feature of the velocity field is that it is monotone with respect to x. Indeed it 

is 



1 + Vy^\ 



since V(x, t) 



^ . 2n _ IEaO-2exp[o-(a; + tM(x,t))] f^^a eicp[a{x + tM{x,t))]\^ ^ ^ 
^""■'^^ ^ ~ E<^exp [a{x + tM{x,t))] ~ \ E<^exp [a{x + tM{x,t))] ) " 

This is known as the entropy condition for the velocity field in the theory of shock waves 
[4][I0]. For M{x,t) it follows that 

d^M{x,t)>0. (19) 

We have seen in [6] that in usual CW model, i.e. with dichotomic spin variables, the spon- 
taneous symmetry breaking associated to the phase transition appears as a shock wave in our 
mechanical analogy. The same happens dealing with our generalized variables. 

Proposition 1. The line {t > tc, 0), with tc = sup ^ is a shock wave for M{x, t), and by putting 
M± = lim^_,o± M{x,t), it is M+ = -M". 

Proof. With a glance to characteristics (fTTl) we notice that x = is a stable point of motioiil]. 
Furthermore for x = all the straight lines do intersect the x-axis in a certain time. Defining 

= sup ^^^°'^^ = supa^M(xo,0), (20) 

XO ^0 XO 

we have that the line {t > t^ 0) is a discontinuity line for M(x, t) since every point on this line is 
an intersection point of characteristics, i.e. it is a shock waves for the velocity field u{x,t). We 
notice from ([20]) that, since we have mf{M'^{x,t)) = 0, it must be 

_ E^a^exp[a{x + tM{x,t))] ^ ^2 
" ~ ^x^ E<,exp [a{x + tM{x,t))] " 

On the other hand, we have that for every time there certainly exists a neighbors of x = where 
the function M(x, t) is smooth. Thus we are allowed to use the Rankine-Hugoniot condition for 
the jump along discontinuity [I][l0] for stating = M^. This last result, coupled with (fT9]l . 
completes the proof. □ □ 



3 Bipartite models 

We are now interested in considering a set of N spin variables, in which is precisely defined a 
partition in two subsets of size respectively Ni and N2. We assume the variable's label of the first 
subset as cxj, i = 1, A^i, while the spins of the second one are introduced by tj, j = 1, ...,N2. 

^It is actually due to the parity of probability measure of a. 
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For each subset all the spins are i.i.d. r.v., with probability measure as discussed above, but in 
principle fj,{a) could be different by /x(t). Of course we have A''i + A''2 = N, and we name the 
relative size of the two subset N2/N1 = ajy- To avoid a trivial behavior of the model, we assume 
that the thermodynamic limit is performed in such a way that a = limiy otv is well defined. 
The spins interact via the Hamiltonian H]y{a,T,hi,h2): 

^ Ni N2 m N2 

i=l j=l i=l j=l 

We notice that spins in each subsystem interact only with spins in the other one, but not among 
themselves. Partition function, pressure and free energy per site for the model are defined 
naturally, in agreement with the previous section: 

AN{P,hi,h2) = -^logZAr(/3, /li,/l2), 
iVi 

/Ar(/3, /ll,/l2) = --Aiv(/3, /li,/l2). 

Remark 4. It should be noticed that, for coherence with already known bipartite models (as i.e. 
the Hopfield model JMi), we choose Ni, instead of N, as the normalization factor inside the free 
energy density and pressure. As we are considering the extensive scaling among the two parties, 
i.e. N2 = a^Ni and limjv^oo «7V = a £ M^, this simply shifts the overall result by a factor 

{i + ayK 

We can also specify the Boltzmann states of our system as 

, , ECTjv,lETjv,5Ar(o",T)exp(-/3i?Ar((T,r,/ii,/i2)) 

As usual, the respective magnetizations of the two systems are 

1 



N2 

N2 



™^ = ^E^^-' (23) 



thus the Hamiltonan reads off as 

H]\f{a, r, h) = —Ni [a]\rm]\rn]\r + hirriN + h2aNn]y] . 

3.1 The free energy in the thermodynamic limit 

In order to reproduce the same scheme of the previous section, let us introduce now the (x, t)- 
dependent interpolating partition function 



ZNix,t)= (24) 
Eo-Er exp A''i f taNTTiNnN + — — —{^"n + ct^n%) + x{mN — aNnj^) + hirriM + /^2aAf?^A^ ) 
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Remark 5. Again we notice that the thermodynamical partition function of the model is recovered 
when t = (3 and rc = 0. 

We can go further and define the action 

ipN{x, t) = ^ log Zn{x, t), (25) 
iVi 

that therefore is just the pressure of the model for a suitable choice of {x,t). Now, computing 
derivatives of ipN{x,t), we notice that, putting Dj\f = ttin — a^riN, it is 

dt<fN{x,t) = -^{D%) {x,t), 
dx^N{x,t) = {DN){x,t), 
dlm{x,t) = ^ (^{D%) - {D^f 

The main difference with respect to the previous case is, instead, the more complicated form 
of the boundary condition, i.e. the action at t = 0. In fact we have that interactions do not 
factorize trivially (in a way independent by the size of the system). It is 

iPn{x, 0) = Alf{P, hi + x) + aNA%{al3, /i2 - x), (26) 

where A]^ is the pressure of the Curie- Weiss model made by A^i a spins, and Aj^ is the same 
referred to the A'2 r spins. Hence, the results of the previous section give us a perfect control on 
the function on the r.h.s. of (f26ll . and we have 



if Nix, 0) = A\l3, hi+x) + aA'^{a(3, h2 - x) + O (^-^^ . (27) 



Thus, again we can build our differential problems for the action ipiy{x,t) 



{ 



dt(pNix,t) + i((9^(^Ar(x,t))2 + ^d^fNix,t) =0 in M X (0, +cx)) 
ipN{x, 0) = A]^{(3, hi+x) + aNAj^{aN(3, h2 - x) on M x {t = 0}, 



(28) 



and for the velocity field DN{x,t) 

{ dtDN{x,t) + DNix,t)dxDN{x,t) + ^dlDN{x,t) = in R x (0, +00) 
\ Dn{x,0) = MN{l3,hi + x) - aNNN{aN(3,h2 - x) on M x {t = 0}, 

whence, as for the boundary condition for the action, we have stated 

Mn{i3, hi+x) - aNNN{l3, /12 - x) = M(/5, hi + x) - aN{aN(3, h2 - x) + O ^ 



(29) 



N 

Remark 6. We stress that our method, due to the existence of the Burger equation for the 
velocity field, introduces by itself the correct order parameter, without imposing it by hands. We 
will back on this point in the last section. 

Remark 7. We have that for each collection of values (/?, a, /ii, /i2), the function Di\i(x,t) is 
bounded^/ N, i.e. the function ipis[{x,t) is Lipschitz continuous. 
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The main difficulty here is that we have a sequence of differential problem with boundary 
conditions dependent by N. Anyway we can replace it with the same sequence of equation but 
with fixed boundary condition, that is the well defined limiting value for — > oo of y?jv and D^. 
To this purpose it is useful the following 



Lemma 2. The two differential problems 

dt^N{x,t) + \{d^ipN{x,t)f + ^dl(pN{x,t) = m M X (0,+cx)) 

(yOjv(a;,0) = A^{l3,hi + x) + aA^{al3,h2 - x) = hN{x) on M x {i = 0}, 



and 



dt(pN{x, t) + \{d^(pN{x, t)f + i^dl^pN{x, t) = in M X (0, +oo) 

(Pn{x, 0) = A^{P, hi+x) + aA^ia/S, h2 - x) = h{x) on M x = 0}, 



(30) 



(31) 



are completely equivalent, i.e. in thermodynamic limit they have the same solution, (pN V and 
(pN — <^ and it is 

\vn-(Pn\ 

Proof. By a Cole-Hopf transform, we can easily write the general form of d]s[{x,t) = \ipN{x,t) — 
(pN{x,t)\ as 



log- 



J:^^dyA{y,{x,t)) 



where we introduced the modified heat kernel A(y, (x, t)) = y ^ exp {^—N [{x — y)'^/2t + h{y)] ) , 

and RN{y) = \h{y) — hN{y)\, with liniM NRn < oo, Vy. Now we notice that because of theorem 
2, it certainly exists an y* such that 

supi?iv(y) = y* and limNR]\r{y*) < oo. 



Hence it is 



y N 



SN{x,t) < l|loge-^^-(«*)| 



= ^[NRNiy*)]<0(^^y (32) 

that completes the proof.D □ 

Of course a similar result holds also for the Burgers' equation for the velocity field Djv- 
So, finally, we must study 



dt(pN{x,t) + \{dx(pN{x,t)f + ^d^(pNix,t) = in R X (0,+oo) 
(Pn{x, 0) = A^{P, ^1 + x) + aA'^{aP, h2 - x) on R x {i = 0}, 



and 



dtDN{x, t) + Dn{x, t)dxDN{x, t) + i^dlD^ix, i) = in M x (0, +oo) 
Dn{x, 0) = M(/?, hi+x)- aN{al3, h2 - x) on R x {i = 0}. 

Now the path is clear, and we can state the following 



(33) 



(34) 
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Theorem 3. The pressure of the generalized bipartite ferromagnet, in the thermodynamic limit, 
is given by: 



A{(3,a,hi,h2 



-a^NM + log E^exp 



a(hi + a(3N 



+ a log E,- exp 



T /12 + /3M 



, (35) 



where, given the well defined magnetization for the generalized CW model respectively for cr and 
T, M{(3,h) and N{(3,h), it is 



Mil3,a,hi,h2) = M{(3,hi- f3M + afiN) 
N{P,a,hi,h2) = N{(3,h2 + PM -a(3N). 



Furthermore it is 



\AN{P,hi,h2)-A{P,a,hi,h2)\ <0[- 



(36) 
(37) 

(38) 



Proof. Theorem 1 gives us the existence and the form of the free solution. We know ^ that the 
free Burger's equation can be solved along the characteristics 



t 

X 



where 



Xq + sZ)(xo,0), 
Z)(xo, 0) = M(/3, hi + Xq) + aN{al3, /i2 - xq), 



(39) 



and it is 

D{x, t) = D{xq{x, t),0) = M{(3, hi+x- tD{xQ, 0)) + aN{a(3, h2 - x + tD{xo,0)). 
Then we can notice that 

Eo-cr exp [a {hi + x + taN)] 



M{P,hi+x-tD{xo,0)) 
N{a(3,h2 - x + tD{xQ,0)) 



exp [a {hi+x + taN)] ' 
ErT exp [r (/i2 - x + tM)] 
Er exp [r ih2-x + tM)] ' 



(40) 
(41) 



which coincide with (l36|) and (|37|l when x = and t = (5. 

At this point we know that the minimum in theorem 1 is taken for y = x — tD{x,t), and, 
bearing in mind the general form of the pressure of CW models, given in the last section, we 
have 



['/5(^>0](x=0,t=/3) 



^D^{x, t) - ^M^i/S, hi+x- tD{xo, 0)) - ^a^N'^{a^, h2 - x + tD{xo, 0)) 



+ log Ea exp [a{hi+x + taN)] + a log E^ exp [r (/i2 - x + tM)] 
= A{(3,a,hi,h2), 



(x=0,t=f3) 



where A{l3,a,hi,h2) is given just by ([35]l . bearing in mind the right definition of M and N. 
Now we must only prove the convergence of the true solution to the free one. But, exactly like 
in theorem 2, equation ([38]) follows by standard techniques, because of the uniform concavity of 



{x - yf 
2t 



+ A^{(), hi+y) + aA'^iP, /i2 - y) 
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with respect to y, assured by theorem 1. In fact we have that, by a Cole-Hopf transform j4], the 
unique bounded solution of the viscous problem is 



2t ' 

and we have, by standard estimates of a Gaussian integral, that 

\(p{x,t) -(pN{x,t)\ <0(^-^ 

i.e. also the ([38]l is proven. □ □ 
Finally, by this last theorem, we can easily write down the free energy of the model: 



1 



a 



f{a,(3,hi,h2) = QiVM- -logE^exp \ cr l^hi+aPNj^ - -logE^exp (^/i2 + /3M 

Remark 8. We stress that when recovering the one party scenario (i.e. a = 0) the model trivially 
reduces to the well known CW in an external magnetic field, with the free energy — /ii) = 
ln2 + lncosh(/?/ii). 

In the last paragraph we will see how expressions like this one can be derived thought a 
minmax principle. 

3.2 The occurrence of a minmax principle for the free energy 

As we have seen in the previous paragraph, the velocity field Diy{x,t) plays the role of order 
parameter for the model. Actually, in perfect analogy with other cases of interest (see for 
instance the last section about generalized ferromagnets, or [6]), the free energy is then obtained 
minimizing (or maximizing, depending on the complexity of the system, i.e. the presence of 
frustration [12]) the action with respect to the order parameter. In bipartite model one has two 
natural order parameters, i.e. each of which referred to the party it belongs to. Prom our study 
of bipartite ferromagnet, we know that the true order parameter is a linear combination of the 
two magnetizations, one for each parties, D = M — aN: What is done by Lax's theorem, for 
example for the free energy, is taking the maximum of on a suitable trial functional [8] 



/(a, /3, /ii, /12) = nmx 



2 2 2 

1 a 

- - log exp [a{hi+(5M + pD)] - - log E^ exp [r (/i2 + a(3N - (3D)] 
fj fj 

This expression is rather unsatisfactory, since not only the order parameter of the model D 
appears, but even the two magnetizations M and A'^. Anyway we can see the model as described 
by two different order parameters, M and N themselves, and in the last expression one should 
take the extremum with respect to both M and A^. Anyway we have that D = M — aN, thus 
maximize D is equivalent to maximize M and minimize A^. We must only rewrite our trial 
functional in terms of M and A^, and we have the minmax principle for the free energy 

r 1 a 

f = min max aMN - - log E<^ exp [a {hi + a(3N)] - - log E^ exp [r (/i2 + /3M)] 
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It naturally arises from the last formula that the free energy is concave with respect to N and 
convex with respect to M, but of course it is uniformly convex along M — aN. Indeed, we have 
that M and N are not independent, but are related by (f36]l and JST]), that is 

^ E^aexp [a{hi + a(3N)] ^ E^rexp [rih^ + pM]] 

exp [a {hi + apN)] ' exp [r (/i2 + (3M)] ' ^ ' 

Remark 9. As for the single party model, we stress that when choosing fi{a) = (1/2) [5(ct + 1) + 
6{a — I)], i.e. dichotomic case, the self- consistent relations reduce to the already known^S^ 

M{l3,hi,a,N) = tanh{hi + paN), (43) 
N{(3,hi,a,N) = tanh(/i2 + /3M). (44) 

However, with respect the model analyzed in IS] it should be noticed that we miss the self- contribute 
inside each equation (i.e. M ^ f{^) o-s well as N ^ fi^) )■ This is ultimately due to the lacking 
of the self-interaction inside each party into the Hamiltonian we are considering. 

These are the true self-consistence relations of the model, analogue to (fTSjl . and we conclude 
that the choice of two different order parameters is redundant, since they are related. One might 
make the choice of putting = N{M) and study the problem using only M as order parameter 
(or viceversa), but, as we have seen it is not so convenient, since a beautiful extremum principle 
does not seem to arise studying the system along the direction of one of the two subsystem^, 
i.e. along M or N. In fact we know, thanks to our technique, that the extremum is taken with 
respect to D. 

Thus actually one has only one degree of freedom, and the minmax principle, although on 
one hand it gives a more satisfactory form of the flow equations, on the other hand only hides a 
more meaningful minimum or maximum principle. This characteristic of bipartite model seems 
to be quite general, and might be extended to other models of interest in future development. 



4 Conclusion 

In this paper we used a mechanical analogy, introduced and developed in [7] [6], for a complete 
resolution of mean field ferromagnetic models with a very general class of spin r.v., i.e. with 
probability measure symmetric and with compact support. The free energy in the thermody- 
namic limit and the phase transition have appeared in our work as, respectively, the solution in 
the limit of vanishing viscosity of a Hamilton-Jacobi equation with diffusion, and the occurrence 
of a shock line for the related velocity field. Moreover, we have applied the same methods to 
the more interesting bipartite systems, made by two different subsystem of spins (a priori of 
different nature), each one interacting with the other, but with no self-interactions. We have 
seen that the thermodynamic limit of the pressure does exist and it is unique and we gave its 
explicit expression in a constructive way. Further, when introducing the Burger's equation for 
the velocity field, our methods automatically "choices" the proper order parameter, which turns 
out to be a linear combination of the magnetizations of the two subsystems, with different signs. 
By this property, we developed an analysis of the minmax principle, pointing out its importance 
relating it to the more classical min/max for the free energy (or, of course, for the pressure) 
for this very simple model. Noticing that the same structure can be recovered for many other 
models of greater interest, like bipartite spin glasses, we plan to report soon about them. 

^This is finally due to the symmetry between the a subsystem and the r one. 
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